Scaling of dynamic spin correlations in BaCu2(Sio.5Geo. 5)207 
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The magnetic dynamic structure factor of the one-dimensional S = 1/2 chain system 
BaCu2(Sio.5Geo. 5)207 is studied in a wide range of energy transfers and temperatures. Contrary 
to previous erroneous reports [T. Masuda et al, Phys. Rev. Lett. 93, 077206 (2004)], the scaling 
properties observed in the range 0.5-25 meV are found to be fully consistent with expectations for a 
Luttinger spin liquid. At higher energies, a breakdown of scaling laws is observed and attributed to 
lattice effects. The results are complementary to those found in literature for other S — 1/2-chain 
compounds, such as KCuFa and Cu-benzoate. 



PACS numbers: 

I. INTRODUCTION 

The two isostructural oxides BaCu2Ge207 and 
BaCu2Si207 are prototypical quasi-one dimensional 
(quasi-lD) quantum S = 1/2 antiferromagnets (AF).^ 
Previous extensive neutron scattering experiments on 
BaCu2Si207 were instrumental in developing an un- 
derstanding of the spin dynamics in weakly-coupled 
5 = 1/2 chainsiSiMi^ More recently, solid solutions of 
type BaCu2(Sii_a;Gea;)207 were recognized as poten- 
tial spin chain systems with random bond strengths^ 
Bulk measurements^*^ and preliminary neutron scatter- 
ing studies^ indicated that the x — 0.5 compound shows 
some very unusual scaling of the static magnetic suscep- 
tibility and the dynamic spin correlation function. It was 
proposed that, due to intrinsic quenched structural dis- 
order, BaCu2(Sio.5Geo. 5)207 is an experimental realiza- 
tion of the much-studied Random Singlet (RS) modeliSiiS 
Unfortunately, when foUowup measurements were per- 
formed on larger and higher quality samples, these con- 
clusions were shown to be erroneous^ In fact, in a 
wide energy range the magnetic dynamic structure fac- 
tor S{q,uj) follows theoretical expectations for disorder- 
free spin chains. In the present paper we report detailed 
temperature-dependent measurements of magnetic exci- 
tations in BaCu2(Sio.5Geo. 5)207 obtained using the new 
large single crystal samples. The focus is on scaling rela- 
tions for the spin correlation function. Convenient mag- 
nitudes of magnetic interactions enable us to cover a wide 
dynamic range huj/nT « 1 — 100. In this, our results 
are complementary to those previously obtained for two 
other 5=1/2 chain systems, namely Cu-benzoate^^ and 
KCuFa (Ref.lH, where ^w/kT « 0.1 - 10. 



BaCu2(Sii_a;Ge2:)207 crystallizes in an orthorhombic 
structure, space group Pnma. The lattice constants lin- 
early depend on Ge content and for x = 0.5 are: a = 
6.917 A, b = 13.28 A and c = 6.944 A. Magnetic proper- 
ties are due to Cu^+ ions that form weakly-coupled anti- 
ferromagnetic chains running along the crystallographic 
c-axis. The in-chain exchange constant is J = 24 meV for 
the silicate {x — 0) and J = 50 meV for the Ge-system 
(x = l)wm The average coupling strength (J) for inter- 
mediate concentrations can be deduced from the posi- 
tion of the Bonner-Fisher susceptibility maximum. It in- 
creases linearly with a;,^'® and for BaCu2(Sio.5Geo. 5)207 
(J) = 37 meV. The corresponding Des Cloizeau-Pearson 
(DCP) zone-boundary energy of magnetic excitations is 
huJzB = Tr{J)/'2 = 58 meV.i'* Nearest -neighbor spin-spin 
separation within the chains is equal to c/2, so the ID AF 
zone center g|| = tt corresponds to / = 1, where {h,k,l) 
denotes a vector in crystallographic reciprocal space. 
Weak interactions between the chains result in long-range 
AF ordering at Tn =9.2 K in the sihcate and Tn =10 K 
in the germanatei However, for intermediate concentra- 
tions, the transition temperature is suppressed^*^ and 
for BaCu2(Sio.5Geo. 5)207, is as low as Tn = 0.7 
Previous theoretical and neutron scattering studies have 
shown that, in weakly coupled spin chains, 3D effects are 
totally negligible at energies fuo > 5KTNrii& so that, bar- 
ring any disorder-induced effects, BaCu2(Sio.5Geo. 5)207 
can be expected to behave like a good ID system for 
hjj > 0.5 meV. Comparing this energy to fiwzB, we see 
that the material presents a conveniently broad energy 
range for testing any predicted ID scaling relations for 
5(q, w) in quantum 5 = 1/2 chains. 
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FIG. 1: Constant-i5 scans measured in BaCu2(Sio.5Geo. 5)207 
using Setup 1 (3-axis spectrometer) at T = 12 K (symbols). 
Heavy solid lines and shaded areas are fits using the Miiller 
ansatz. Dashed lines are Lorentzian fits. All fitting functions 
are convoluted with the known experimental resolution. 



II. EXPERIMENTAL 

In the present study we used 6 BaCu2(Sio.5Geo. 5)207 
single crystals of total mass 15 g, grown using the float- 
ing zone method. Most crystals were twinned in such a 
way that the a and c axes of the two crystallographic 
domains were interchanged. Such a growth habit is con- 
sistent with awe, but complicates the measurements, 
as it involves the chain axis. Individual crystals were 
co-aligned to form a compound sample with a total ef- 
fective mosaic of 1.7°. Inelastic neutron scattering data 
were collected in two separate series of experiments. Ex- 
perimental Setup 1 employed the HBl thermal 3-axis 
spectrometer installed at the High Flux isotope Reactor 
at ORNL. Final neutron energy was fixed at 14.7 meV. 
A Pyrolitic Graphite (PG) monochromator was used in 
conjunction with a PG filter installed after the sample. 
The colfimation setup was 48' - 40' - 40' - 240'. The 
sample was mounted with the (a, c) scattering plane hor- 
izontal, making crystallographic wave vectors (/i, 0, 1) ac- 
cessible for measurements. Sample environment was a 
standard He-fiow cryostat. All data were corrected for 
the wavelength-dependent efficiency of the beam moni- 
tor (due to higher-order beam contamination) and the 
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FIG. 2; (Color online) Time of flight spectrum collected in 
BaCu2(Sio.5Geo. 5)207 at T = 5 K using Setup 2, in projec- 
tion onto the c (chain) axis. Solid and dashed white lines 
are ID antiferromagnetic zone-centers for the two crystallo- 
graphic domains present in the sample, respectively. 



vertical focal length of the fixed-curvature monochroma- 
tor. It was, in part, an incorrect application of these 
corrections, that led to erroneous conclusions in Ref. 0. 
The second series of measurements were performed on 
the MAPS time-of-flight (TOF) chopper spectrometer in- 
stalled at the ISIS spallation neutron facility at RAL. In- 
cident energy was fixed at either 50 meV (Setup 2) or 
100 meV (Setup 3). The (a,c) plane of the crystal was 
mounted parallel to the preferential scattering plane of 
the instrument. All TOF data were integrated along the 
b* crystallographic direction in the range — 1 < fc < 1. 
Sample temperature was controlled between 5.5 K and 
300 K using a closed-cycle refrigerator. 



III. RESULTS AND DATA ANALYSIS 
A. 3-axis data 

In 3-axis experiments the data were collected in 
constant-i? scans. Due to the presence of two sets of 
crystallographic domains, any scattering at wave vector 
(/i, fc, I) is observed simultaneously with that at (Z, fc, h). 
A steep dispersion of magnetic excitations along the 
chains implies that, at low energy transfers, all relevant 
scattering is concentrated near the ID AF zone-centers l- 
odd. For this reason, scans performed at Z = effectively 
pick up the signal from only one of the two domains. 
While only half of the sample contributes to scattering, 
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FIG. 3; (Color online) Time of flight spectrum collected in BaCu2(Sio.5Geo. 5)207 at T = 5 K using Setup 2, in projection onto 
the c (chain) axis. Solid and dashed white lines are ID antiferromagnetic zone-centers for the two crystallographic domains 
present in the sample, respectively. 



this geometry is optimal for achieving a high wave vector 
resolution along the chains. A typical constant-i? scan 
of this type measured using Setup 1 at huj = 12 meV 
is shown in Fig. An alternative approach is to re- 
strict data collection to q = {h ± x,k,l ± x), where x 
is arbitrary and h and I are integers. At these positions 
the scattering in both crystallographic domains occurs 
at wave vectors that are equivalent in the ID sense. The 
presence of two domains still has to be carefully taken 
into account when calculating the spectrometer resolu- 
tion function. This second approach is the only option 
of high energy transfers, where in each domain the scat- 
tering is spread out over a wide g-range, and necessarily 
interferes with that from the other domain. Constant-i? 
scans collected in this mode are shown in Fig. ^,c. A 
constant background was subtracted from all 3-axis scans 
shown. 

At low temperatures the magnetic dynamic struc- 
ture factor of disorder-free quantum S — 1/2 
chains can be accurately approximated by the Miiller 
AnsataS^aiLi&i2i2a*2i (MA). For BaCu2(Sio.5Geo.5)207, 
m Ref.i it was claimed that in constant-i? scans the mag- 
netic scattering appears broader than the MA expecta- 
tion, and is more consistent with Lorentzian line shapes. 
It was since realized that the elongated Lorentzian-type 
"tails" are a result of background contamination.^-'^ An- 
alyzing the new data we find that the MA peak profiles 



actually work for BaCu2(Sio.5Geo. 5)207 rather well. For 
each measured constant- i? scan the MA dynamic struc- 
ture factor was numerically folded with the known spec- 
trometer resolution function and fit to the experimen- 
tal data. The DCP spin wave velocity v was fixed at 
7r(J)c/4 = 404 meVA. The only variable parameter was 
an overall scaling factor. Excellent fits were obtained 
at all energies (solid lines and shaded areas in Fig. 
Simulations using Lorentzians convoluted with the ex- 
perimental resolution function are shown in dashed lines. 
Lorentzians provide satisfactory fits only at low energy 
transfers, where the finer spectral features are masked 
by the limited experimental resolution anyway. At high 
energies deviations from Lorentzian line shapes become 
obvious, and a flattening at the intensity maximum re- 
sulting in typical "top-hat" shapes is clearly observed 
(Fig.Ct). 



B. TOF data 

Specifics of the TOF geometry make separating the 
signals originating from the two crystallographic domains 
less straightforward. Typical data collected using Setups 
2 and 3 are shown in Figs. [3 and |3| in projection onto 
the c-axis of one of the domains (domain A). Vertical 
streaks of scattering at the ID AF zone-centers I — ±1 for 




FIG. 4: Typical constant-i? cuts of the TOF data for 
BaCu2(Sio. 5060.5)207 collected using Setup 2 at T = 4 K. 
Solid lines are Gaussian fits. A flat background has been sub- 
tracted. 



FIG. 5: Temperature evolution of a constant-iJ cut measured 
for BaCu2(Sio.5Geo. 5)207 using Setup 3. Solid lines are Gaus- 
sian fits. A flat background has been subtracted. The shaded 
area is due to scattering by phonons. 



this domain are clearly visible and highlighted by vertical 
solid lines. However, they are intersected by the ID AF 
zone-centers of the 2nd domain (domain B), represented 
by the dashed curves. The corresponding magnetic scat- 
tering from domain B appears as a series of arcs. In order 
to separated it from domain-A scattering one has to rely 
on certain assumptions and approximations. At small 
energy transfers, where all the scattering is concentrated 
in the vicinity of the ID zone-centers, one can simply re- 
gard all data collected sufficiently far from domain B ID 
zone-centers as belonging to domain A. Constant-energy 
cuts obtained from Setup 2 data using this assumption 
are shown in Fig. ^ (constant background subtracted). 
Compared to 3-axis data, the experimental wave vector 
resolution of the TOF setup is insufhcient to resolve any 
fine structure in the scattering profile. 

An alternative approach is to select a geometry in 
which domain B ID zone-centers cross that of domain 
A at an almost constant energy, as in the case of Z = 1 
in Fig. O For a constant-i? cut extracted from such a 
data set, the magnetic scattering by domain B will, to 
a good approximation, correspond to a constant Z-value. 
Since ID magnetic scattering is independent of h and fc, 
it will amount to a constant background contribution. 
Constant-energy scans extracted from Setup 3 data us- 
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FIG. 6: Scaling of S[uj) measured in BaCu2(Sio. 5060.5)207 
(symbols). Solid circles: Setup 1 (3-axis), T = 12 K. Open 
circles: setup 2 (TOF), T = 5.5 K. Squares: Ref0, re-analized 
thermal-neutron 3-axis data, T = 12 K. Triangles: Ref0, 
cold-neutron 3-axis data, T — 1.5 K. Solid line: Luttinger 
spin liquid at T — > 0. The dashed curves are guides for the 
eye. Inset: same plot for hujS{ijj). 
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over the entire g-range, but, on a lattice, it should be 
limited to the vicinity of the ID AF zone-center. Taking 
the integral entirely eliminates the effect of wave vector 
resolution. In our case of BaCu2(Sio.5Geo. 5)207, S{ijj) 
was determined from the cross section fits to individual 
3-axis scans q-scans (solid lines in Fig.^. For the TOF 
data we employed empirical Gaussian fits to constant-£' 
cuts, as shown in solid lines in Figs. 01 and [S] The bulk 
of the low-temperature data for S{iu) are shown in Fig.|B| 
(circles). Re-analyzed data from Ref. are plotted as 
solid squares (thermal neutron measurements) and tri- 
angles (cold neutrons). The inset shows the same data 
scaled by the energy transfer. All S{uj) data collected in 
this work at different temperatures, after an appropriate 
scaling (see below) , are plotted in Fig. [7| Solid symbols 
are for energy transfers fuv < 25 meV, while open sym- 
bols represent higher energies. 



FIG. 7: Measured temperature scaling of S(tj)/(n(j -I- 1) for 
BaCu2(Sio.5Geo. 5)207. Open symbols correspond to energy 
transfers above 25 meV. The solid line is the approximate 
scaling function for a Luttinger spin liquid^ Dashed lines are 
guides for the eye and emphasize the breakdown of scaling at 
high energies. 



ing this procedure are shown in Fig. [S| A temperature- 
induced broadening of the scattering profiles is clearly 
visible. Data analysis at high temperatures is compli- 
cated by the increasing phonon contribution to the back- 
ground (Fig. and shaded areas in Fig. [Sj. For this 
reason, the measurements were limited to T < 100 K. 



C. Scaling quantities 

The central purpose of this study was investigating the 
energy and temperature scaling of the spin correlations 
in BaCu2(Sio.5Geo. 5)207. Much of the existing experi- 
mental work concentrated on the scaling of the dynamic 
structure measured at the ID zone-center S{7r,uj)m^^ 
This quantity is defined exclusively by long wavelength 
AF spin correlations. In this regime quantum spin chains 
are equivalent to a Luttinger spin liquid. The exist- 
ing exact analytical predictions for this model^^S make 
a comparison with theory straightforward. The problem, 
however, is that measuring S{tt,uj) requires wave vector 
selectivity: the experimental wave vector resolution 6q 
must be small compared to the g-width of scattering at 
a particular energy transfer. The latter can be estimated 
as huj/v, where v is the DCP spin wave velocity. Thus, 
there will always be a lower limit on energy transfer for 
the measurements, determined by the experiemntal res- 
olution: hujmin ^ Sqv. For our 3-axis setup, the practical 
lower limit is about 20 meV, and is still higher for the 
TOF experiments. A more robust approach is to measure 
the g-integrated structure factor S'o(w) = / S{q,uj)duj^- 
In the hydrodynamic limit the integral is to be taken 



IV. DISCUSSION 

The finite-temperature scaling for the dynamic struc- 
ture factor of a Luttinger spin liquid is given byi^ 



^/ - N / N 1 f huj vq 

S{q,uj) a (n^ + l)_Ff — 



(1) 



Fix V) - Im M U 47r ) (2) 

Fi-^y) ^"l^r(|-^^)r(|-z^);' 

where {n^ + 1) = (l - e-'""^/^'^^)) "\ For an 5 = 1/2 
Heisenberg spin chain, this scaling applies to AF spin 
correlations near a ID AF zone-center q — qo, so that 
q = q—qo- The integral over q needed to derive the scaling 
for S{lu) cannot be easily taken analytically. However, it 
can be very accurately approximated bytii 

huj 



S{uj) ^ tanh 



2kT 



(3) 



To a very good approximation, then, S{uj) — const at 
hui ^ kT. The same result can be obtained by integrat- 
ing the MuUer Ansatz, which at small q and Hio <^ huzB 
coincides with the T — > asymptotic form for the scaling 
function 121 

The low-temperature data plotted in Fig. show that 
in BaCu2(Sio.5Geo. 5)207, S{lj) is indeed practically en- 
ergy independent up to about 25 meV energy transfer 
(solid lines), i. e., up to about half of the zone-boundary 
energy fkJzB ■ In this energy range the measured temper- 
ature dependence of S{uj) is consistent with the scaling 
form|31 as shown in Fig.Q solid line. This plot, with a dy- 
namic range for Hlo/kT from 1 to 100, is complimentary 
to that shown in Ref. 17 for Cu-benzoate with Hlo/kT in 
the range of about 0.1 to 10. Above 25 meV, the mea- 
sured S{uj) starts to decrease and ujS{lo) deviates from 
linear behavior, even at low temperatures (Fig.|2| dashed 
lines). This breakdown of scaling is observed at all tem- 
peratures, as shown by the open symbols and dashed lines 
in Fig. 13 
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The observed breakdown of scaling for S{uj) at about 
half the zone-boundary energy is rather interesting, if not 
unexpected. Work on KCuFa has shown that for 5(7r, uj) 
the scaling relation ^ holds at much higher energies, up 
to about 140% of Iiujzb-^'^ Of course, the difference is due 
to the fact that Eq. ^ is derived for the continuum limit 
\q\ ^ 1. Unhke S'(7r, w), S{lo) includes scattering at wave 
vectors —Hlu/v < q < —hui/v, and thus rapidly exceeds 
the domain of validity of Eq. ^as hiu/v increases. Note 
that the existing S{u!) scaling data for Cu-benzoate also 
extend to about half of Hluzb- 



pounds are useful S = 1/2 chain materials, with a con- 
tinuously adjustable exchange constant, is fully justified. 
On the other hand, the suppression of the ordering tem- 
perature, the anomalous low-temperature increase in the 
bulk susceptibility,® recent ESR data^'^, and preliminary 
neutron spin echo measurements^^ clearly show that, at 
lower energies (longer time scales), quenched disorder 
does become relevant. Understanding the scaling of spin 
correlations in this regime is a challenge for future stud- 
ies. 



V. CONCLUDING REMARKS 
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